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Abstract 

The presence of random fields is well known to destroy ferromagnetic order in Ising systems in two dimensions. When 
the system is placed in a sufficiently strong external field, however, the size of clusters of like spins diverges. There is 
evidence that this percolation transition is in the universality class of standard site percolation. It has been claimed that, 
for small disorder, a similar percolation phenomenon also occurs in zero external field. Using exact algorithms, we study 
ground states of large samples and find little evidence for a transition at zero external field. Nevertheless, for sufficiently 
small random field strengths, there is an extended region of the phase diagram, where finite samples are indistinguishable 
from a critical percolating system. In this regime we examine ground-state domain walls, finding strong evidence that 
they are conformally invariant and satisfy Schramm-Loewner evolution [SLE^) with parameter k — Q. These results 
add support to the hope that at least some aspects of systems with quenched disorder might be ultimately studied with 
1 the techniques of SLE and conformal field theory. 
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The random field Ising model (RFIM) is one of the 
earliest studied and simplest disordered systems showing 
non-trivial and glassy behavior pQ [5] . It has a number of 
important realizations in nature, including diluted antifer- 
romagnets in a field and binary liquids in porous media 
[2]. Through its long history, researchers have managed 
to gain a reasonable understanding of the critical behav- 
ior, although this progress has been neither straight nor 
smooth, and many questions remain unanswered [5]. It 
is known, for example, that the RFIM in two dimensions 
(2D) lacks ferromagnetic order OE]. Even at zero tem- 
perature it remains in the paramagnetic state for non-zero 
disorder. Numerical ground-state calculations have shown, 
however, that even in the absence of a thermodynamic 
transition there exists a geometric transition at which the 
size of the spin clusters diverges in a manner bearing many 
similarities to classical site percolation [5l [6] . While this 
transition is rather clearly established in the presence of 
an external field, it has been argued that a similar per- 
colation phenomenon can also be observed in the absence 
of an external field for sufficiently small disorder [6l [?]■ 
Here, we re-investigate the zero-field behavior with large 
scale ground-state calculations, focusing on the possible 
percolation phenomenon. 

The observed relations to classical site percolation at 
this (non-zero or zero field) geometrical transition moti- 
vate further questions of how far the similarities go. In- 
terfaces in two-dimensional percolation satisfy Schramm- 
Loewner evolution (SLE) [Slinj: but this property relies on 



conformal invariance which is not conserved in the pres- 
ence of disorder. Recently, however, there have been sug- 
gestions that the domain walls of certain other disordered 
systems satisfy SLE. This intriguing possibility implies 
that conformal invariance, broken by disorder, is restored, 
at least at criticality where relevant length scales diverge. 

Schramm-Loewner evolution is a method for construct- 
ing a statistical ensemble of curves in the plane from one- 
dimensional Brownian motion, thus classifying curves with 
only one parameter, the diffusion constant k [TU]. Char- 
acteristic interfaces in many physical systems have been 
shown (in some cases rigorously) to satisfy SLE^. These 
include percolation (k — 6), self avoiding walks {n — 4/3), 
as well as spin cluster boundaries {k = 3) and Fortuin- 
Kasteleyn cluster boundaries (k = 16/3) in the Ising model. 
A number of numerical studies have found interfaces in cer- 
tain disordered systems consistent with SLE, in particular 
the 2D Ising spin glass [TTl FT^ . the Potts model on dy- 
namical triangulations [13] , the random bond Potts model 
[13], and the disordered solid-on-solid model 15 . Here we 
extend this list to include the 2D random field Ising model. 

In Sec.[l]we introduce the RFIM and Schramm-Loewner 
evolution in more detail, and discuss the approach used 
here for determining ground states of large samples. Sec- 
tion [2] is devoted to an investigation of the critical behav- 
ior of the RFIM near the geometric transition, focusing 
on the behavior in zero external field for small disorder. 
In Sec. |3j we report the results of our tests of the corre- 
spondence between interfaces in the RFIM and Brownian 
motion implied by SLE. 
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1. The model and method 

We consider the random field Ising model in two di- 
mensions with Hamiltonian 



T~L — —J / ^ SiSj — y ^ hiSi 



(1) 



Here, the spins Si = ±1 are located on the sites of a square 
lattice and interact ferromagnetically with nearest neigh- 
bors. The local fields hi are quenched random variables 
drawn from a normal distribution with mean H and stan- 
dard deviation A. Since, at zero temperature, only the 
ratio J/ A is relevant, we take J = 1 for simplicity. The 
spin-spin interaction J induces a correlation between the 
spins resulting in spin clusters which are compact up to a 
length ^j). Above this scale the clusters are fractal objects, 
the magnetization is zero and the system is paramagnetic. 
As the randomness A is decreased, the breakup length ^(, 
increases. At and above three dimensions ^f, diverges at the 
thermodynamic phase transition, below which the system 
is ferromagnetic. In two dimensions no thermodynamic 
phase transition exists, and ^b diverges only at A = 0. It 
has been argued that, for iJ = 0, the breakup length scale 
^f, increases with decreasing A as [3] 



.A/A- 



(2) 



Though there is no thermodynamic transition in 2D, the 
linear extent of the largest clusters diverges for sufficiently 
large H. In most aspects, this divergence appears to be 
consistent with standard site percolation [H O [7] . It has 
been suggested that the divergence occurs even dX H — Q 
if A is below a critical value 1^ . It is the characteristics 
of clusters of aligned spins and their boundaries at this 
geometric transition which we focus on in this study. 

We restrict our investigation to ground-state spin con- 
figurations, which can be efficiently constructed through 
a mapping to the well known minimum cut (or maximum 
flow) problem in graph theory [TBI [17]. Consider a directed 
graph with A^ + 2 vertices, and edges (i, j) furnished with 
weights Cij. The minimum (s, t) cut is given by a subset of 
the edges of the graph, whose removal disconnects vertices 
s and t, such that the sum of the weights of the cut edges 
is minimal. Using the variables Xi , which are 1 if vertex 
i is connected to s and otherwise, the total weight of a 
cut can be represented as 



N 

C({x}) = Y. ^^(1 



•^j ) ^i 



(3) 



With the identification of Xi with spin variables (except- 
ing the vertices s and t) , and an appropriate choice of the 
weights Cij, this function can be made to precisely match 
the RFIM Hamiltonian (jTl). The minimum cut separating 
the graph into vertices connected to s and vertices con- 
nected to t gives the minimum-energy way of cutting the 
RFIM lattice into clusters of up and down spins, and thus 



corresponds to the ground state of equation [T] The choice 
for the edge weights is, for i ^ (s, t) 



0, 
4J, 



i >j 
else 



(4) 



For the edges connecting the "spin" vertices to s and t, 
the result is given in terms of the quantity Ui — —hi — 



5E,(c, 



i) as 



0, 



Ui> 
else 



Cit 



Ui > 



0, else 



(5) 



Cis and Cti are taken to be zero for all i. Here, we use a 
fast algorithm for solving the minimum cut problem based 
on the idea of "augmenting paths" fTS", 19] . The worst case 
scenario for the running time of this algorithm is an unim- 
pressive 0{N'^) (or more generally V^E where V is the 
number of vertices and E is the number of edges), how- 
ever, the algorithm was designed to optimize the typical 
case. The optimization was carried out for vision and im- 
age analysis problems, but even for the graph structure 
corresponding to the RFIM ground state calculation the 
running time is proportional to A^ for the samples consid- 
ered here. In practice, the maximum system size is limited 
more by computer memory constraints than by time. 

Schramm-Loewner evolutions are defined in terms of a 
family gt of conformal maps which take, formally, the up- 
per half plane minus the curve jt (parametrized by "time" 
t) to the upper half plane, gt : H \ 74 — > H. This map can 
be defined (using complex notation and suitable normal- 
ization) in terms of the differential equation 



dt 



9t{z) -^t 



(6) 



where ^t is the unique driving function for the curve 74. 
If the curve (or the process generating the curve) satis- 
fies SLEk, then ^j will be a Brownian motion with zero 
mean and variance nt. Numerically, rather than solving 
the differential equation, the map gt is instead pictured 
as a series of maps gi which iteratively remove a small 
section from the beginning of the curve. To calculate the 
driving function from a given curve, the incremental map 
gi is approximated using a vertical slit map |20j 



g,{z) ^ i^/^{^^'^^)^^lAti + ^, 



(7) 



The parameters ^i and Ati are determined from the coor- 
dinates of the curve segment to be removed through the 
relations £,i = Xi^i-i and Ati = yf,i_i/4. More specifically, 
Xi^i-i and yt.i-i are the coordinates of the i'th segment 
of the curve after undergoing the i — \ successive maps 
gt-i o . . .0 gi. The parameter S,i is the value of the driving 
function ^t sampled at time ti — X^iXi ^^j- The complex 
square root in equation [7] is calculated, as usual, with the 
branch cut along the negative real axis. 

This is an iterative process in which the coordinates of 
the interface are successively updated for each step, and 
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Figure 1: Schematic representation of the dependence of the location 
He of the percolation transition on the random-field strength A. For 
H < He the spin clusters have finite extent. For H > He there exists 
at least one infinite cluster. 



thus the computational complexity is 0{Lj) where Lj is 
the length of the interface. If L/ grows like ^ A^'^//^ (with 
df = 7/4 expected for percolation) in the number of spins 
N, the resulting computational complexity is 0(7V^/^), 
which is significantly slower than the ground-state calcu- 
lation which is 0{N) on average. We use a fast implemen- 
tation of this "zipper" algorithm [5T], in which blocks of 
multiple slit maps are approximated by a Laurent series. 
Treating blocked maps in one step dramatically speeds up 
the calculation such that it scales on average as Lj-^. The 
loss in accuracy from the approximation is minimal. 

If the driving function ^t can be shown to be Brow- 
nian motion then the curves satisfy SLE^. In practice, 
the finite size of the lattice and the zipper algorithm in- 
troduces correlations between the increments of ^i and in 
their associated time values ti. In particular the distribu- 
tion of time steps ti — ti_i is highly non-trivial and has 
significant correlations. However, it should be emphasized 
that correlations in the times ti at which S.t is sampled 
does not imply correlations in the underlying continuous 
driving function ^f 

2. Phase diagram and behavior at zero field 

For very large disorder, A 3> Jz, where z is the lat- 
tice coordination number, the interactions between the 
spins play no role and the system reduces, trivially, to the 
classical site percolation model. Each spin is determined 
solely by the independent random variable hi. Identify- 
ing spin up with "site occupation", the site occupation 
probability p is simply the probability for hi to be pos- 
itive: p = ^erfc{— i//(-\/2A)}. For smaller disorder the 
spin-spin interaction of the RFIM complicates the anal- 
ogy with percolation. However, it has been demonstrated 
that there exists a line of critical external fields Hc{A), as 
pictured schematically in Fig. [TJ for which observables like 
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Figure 2: Crossing probabilities for spin clusters of the 2D RFIM 
on rectangular domains of aspect ratio e~^'^ for a number of system 
sizes A'^ = e~^'^L'^ at H = 0. Top panel: probability of crossing of 
clusters of up spins in the vertical direction. For A decreasing from 
infinity, the crossing probability approaches the exact percolation 
crossing value |22j . indicated by the horizontal line. At very small 
disorder, the breakup length ^f, is larger than the system size and 
the systems are effectively ferromagnetic. Inset: enlarged view of the 
plateau region for the largest system size. The bottom panel shows 
1 — TTtot, the total crossing probability. If there was a transition at 
H = 0, the lines would have to cross around the transition point. 



the crossing probability and the fractal dimension of the 
spin clusters maintain the percolation values [6j. 

For very large disorder the line of critical external fields 
is found from the critical site percolation probability pc 
{pc ~ 0.5927 for the square lattice) 



i?c(A > z) « -AV2erfc"^{2pJ. 



(8) 



For small disorder the behavior is less well understood. 
7Jc(A) decreases as A decreases, approaching the limiting 
value -ffc(A = 0) = 0. It has been claimed Oil] that at a 
finite disorder strength, A^ « 1.65(5), the curve becomes 
identically zero i?c(A < A^) = 0, cf. Fig. [T] We tested 
this claim by looking at the behavior of the system vary- 
ing A at iJ = 0. We calculated the probability that a 
connected cluster of up spins crosses a rectangle of aspect 
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Figure 3: Distribution of the random walk at several different 
"times" t as extracted from the RFIM interfaces at A = 1.65. The 
lines are normal distributions with zero mean and variance Kt with 
K = 6. 



ratio e~^/^ in the vertical (shorter) direction (this aspect 
ratio is chosen to be different from 1/2, but the specific 
value is somewhat arbitrary) . In the upper panel of Fig. ^ 
we present the results of these calculations. The crossing 
probability n for large A is zero (as expected for a square 
lattice with pc > 0.5). As the disorder is decreased, tt 
approaches the exact percolation value, which is indicated 
by the horizontal line in the plot. At very small disorder, 
when the breakup length scale ^b becomes comparable to 
the system size L, the system appears ferromagnetic and 
the crossing probability falls away from the plateau value. 
To separate the approach to the plateau from the fi- 
nite size effects we also considered another quantity, the 
probability iTtot that there exists a connected cluster of 
either up or down spins crossing either horizontally or 
vertically. The small disorder limit of this quantity is 
7rtot(A — )• 0) = 1 for both the onset of percolation and the 
ultimate ferromagnetic (and finite size) ordering. If there 
is a percolation transition at finite disorder (Ac > 0), in 
the thermodynamic limit L — oo the curve -Ktoti^) will 
be a step function. Hence, the curves for finite systems 
approaching this step function will intersect at Ac. As is 
apparent from our data for 'Ktot shown in the lower panel 
of Fig. [2] such a crossing is not observed at least down to 
A « 1.45, significantly below the previously conjectured 
value of Ac ^ 1.65. The study of even smaller disorder 
strengths, while ensuring L ^ ^f,, is preempted by the ex- 
ponential growth of the breakup length ^(,. Although it 
seems very unlikely that the system undergoes a true per- 
colation phase transition for H = 0, there exists a large 
plateau region in which the system appears to be at crit- 
ical percolation at iJ = 0, even up to very large system 
sizes. Presumably, this holds true even for system sizes 
occurring in experimental realizations of the RFIM. 



3. Schramm-Loewner evolution 

Finally, we tested directly the conformal mapping upon 
which SLE is based. We studied the statistics of 10 000 
interfaces generated in a half disc. This geometry is used 
to optimally mimic the full half plane. The interface is 
initiated at the origin (the center of the flat edge) by two 
fixed spins and is considered ended when it touches the 
curved boundary. We find that the variance of the driving 
function calculated from the interfaces using the method 
described in Sec. [l] is k = ((6 - (6))^) A = 6.09 ± 0.09, 
and the normalized mean is ^ = {£,t) /Vki = 0.01 ± 0.01. 
The agreement with SLE is good, but is only expected to 
be perfect at criticality. The difference from the expected 
percolation value k = 6 which is perhaps just visible here is 
due to the calculations being carried out at A = 1.65. The 
position distribution of the resulting stochastic process at 
several fixed times is shown in Fig. [3] along with curves 
representing the expectations for a perfect random walk. 

We have focused here on the case of zero external field, 
where we have shown, using significantly larger system 
sizes than had been accessible before, that there is no 
percolation phase transition for iJ = 0, at least down to 
A = 1.45. Due to the fact, however, that the percolation 
transition line comes exponentially close to H = Q for small 
random-field strengths A (cf. Fig.jT]), there exists a plateau 
region where even at -ff = the behavior appears nearly 
indistinguishable from criticality. In this regime we have 
tested the applicability of Schramm-Loewner evolution to 
the RFIM and found good agreement. Further studies at 
non-zero external field and using an array of further tests 
for consistency with SLE confirm this result [23] . 
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